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INTRODUCTION 


One of the design requirements of the Space Transportation System (STS) vehicles 
dictated that the vehicles be capable of controlled fliqht during entry through the 
entire flow regime from f ree-molecule through hypersonic to subsonic flow. The re- 
sulting vehicle resembles in many ways a conventional aircraft in that it is a winged 
spacecraft with elevons, vertical tail, rudder, and a body flap trim device. The 
elevons are used both for longitudinal pitch control, much like elevators, and for 
lateral control, like ailerons. These aerodynamic control surfaces are augmented 
with onboard reaction control pitch and yaw jets which are necessary for the low 
dynamic pressure regime. 

Large quantities of wind-tunnel data were gathered during the design of the 
space shuttle. The accumulated data base describes the assumed aerodynamic charac- 
teristics of the shuttle over a wide range of flight conditions. This data base, 
published in reference 1, will be called herein the preflight or data book values. 

Nine shuttle flights (STS-6, 7, 8, 11, 13, 17, 24, 26, and 30) were flown by the 
shuttle vehicle Challenger, Since no additional flights of this vehicle are pos- 
sible, the purpose of this paper is to summarize the extraction of lateral stability 
and control derivatives from lateral maneuver data obtained during entry of the 
Challenger into the atmosphere. The results presented herein constitute part of the 
research conducted at Langley Research Center to analyze the aerodynamics of the 
shuttle vehicle (refs. 2-9). 

Lateral maneuver data were available for six of the nine fliqhts. During two of 
the fliqhts (STS-24 and STS-30) no data were measured; for STS-26, the data were mea- 
sured but not available for analysis. Of the remaining six flights, 33 lateral ma- 
neuvers specifically designed for parameter extraction (called a Programmed Test 
Input or PTI) were performed on five flights; on the sixth flight (STS-17), five 
other lateral maneuvers were analyzed. The 38 lateral maneuvers constitute the data 
base for the present study. Because of safety constraints, the maneuvers are not 
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optimal for parameter extraction; however they are the best available fliqht data for 
the purpose of this study. The fliqht extracted values are compared to the prefliqht 
values of reference 1. 

SYMBOLS 

a y acceleration in y-direction, q units 

b wing span, m 

rolling-moment coefficient, ^x/^w^w 

C 0 .c aerodynamic moments for trimmed fliqht 

x , o ri , o 

C n yawing-moment coef f icients , 

C v _ aerodynamic force for trimmed fliqht 

1,0 

Cy lateral-force coefficient, F y/^ s w 

e vector of measurement error 

F vector function representing equations of motion 

g acceleration due to gravity, 9.81 m/sec 

G vector function representing measurement equations 

I x ,Iy,I^,I x ^ moments of inertia 
J cost function 

k number of data points 

L likelihood function 

m mass, kg 

p roll rate, rad/sec 

q pitch rate, rad/sec 

2 

dynamic pressure, pv /2, Pa 

Q vector of unknown parameters 

r yaw rate, rad/sec 

R measurement noise covariance matrix 

S wing area, m^ 
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t 


time, sec 


a velocity along X-body axis, m/sec 

U input vector 

v velocity along Y-body axis, m/sec 

V airspeed, m/sec 

w velocity along Z-bodv axis, m/sec 

X vector of states 

X , Y , Z longitudinal, lateral, and vertical body axes 

Y vector of outputs 

a angle of attack, rad 

6 sideslip angle, rad 

6a aileron deflection, rad 

5r rudder deflection, rad 

5RCS RCS control term, number of jets firing 

0 pitch angle, rad 
roll angle, rad 

4> bias on roll rate, rad/sec 

o 

Subscripts: 

1 quantity at ith time 

M measured quantity 

p,r rotary derivatives 

0 static derivatives with respect to 0 

6a,5r,6RCS control derivatives with respect to indicated quantity 
t trimmed value 

Matrix exponents: 

T transpose of matrix 

-1 inverse of matrix 
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Mathematical notation: 


A estimated quantity when over symbol 

derivative with respect to time when over symbol 
V gradient operator 

Abbreviations : 

ACIP Aerodynamic Coefficient Identification Packaqe 

BET Best Estimated Trajectory 

IMIJ Inertial Measurement Unit 

DFI Development Flight Instrumentation 

LaRC Langley Research Center 

MMLE3 Modified Maximum Likelihood 

PTI Programmed Test Input 

RCS Reaction Control System 

RGA,AA Rate Gyro Assembly, Accelerometer Assembly 

STS Space Transportation System 


Test Vehicle 

The orbiter configuration is shown in figure 1 and key physical characteristics 
are given in table 1. The thick, double delta wing is configured with full span ele- 
vons, comprised of two panels per side. Each elevon panel is independently actuated. 
All four panels are deflected symmetrically as an elevator for pitch control, and 
left and right elevons are deflected differentially as an aileron (6a) for roll 
control. 

The body flap is used as the primary longitudinal trim device. The elevons are 
programmed in conjunction with the body flap to follow a set schedule to provide the 
desired aileron effectiveness. 

The vertical tail consists of the fin and a split rudder. The rudder panels are 
deflected symmetrically for yaw control and are separated to act as a speed brake to 
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provide for subsonic energy modulation. The speed brake opens fully (87.2 degrees) 
just below Mach 10 and then follows a predetermined schedule until Mach 0.9 is 
reached. The rudder is not activated until Mach 5. 

Stability augmentation is provided by the aft reaction control system (RCS) 
jets, with the forward jets reserved for on-orbit attitude control and for aborts. 

The aft yaw jets are active until Mach 1, while the pitch and roll jets are termi- 
nated at a pressure of 20 and 10 pounds per square foot, respectively. Additional 
details of the shuttle vehicle and its systems are given in reference 1. 

Maneuvers 

During flights STS-6, 7, 8, 11, and 13, especially designed Programmed Test 
Input (PTI) maneuvers were performed to obtain data for use in extracting aerodynamic 
parameters. These maneuvers were performed to obtain data at specific points during 
the descent trajectory. The test points were chosen so that aerodynamic parameters 
could be determined along the descent trajectory to verify the aerodynamic model ob- 
tained from the wind tunnel tests. This verification procedure adds confidence to 
the assumed aerodynamics of the shuttle where there is agreement and points to areas 
of potential inaccuracy where there is no agreement. 

The actual forms of the inputs to be performed were developed using a shuttle 
simulation to generate responses for various inputs and then extracting parameters 
from these responses. The control inputs that gave the best definition of the param- 
eters of interest were then used for the flight tests. In spite of the care taken to 
design effective inputs and because the automatic control system was active, the con- 
trols were coupled and the resulting responses were reduced in magnitude and corre- 
lated with each other and the control inputs. This led to identif iability problems 
and correlation of parameters during the extraction process. Additional details on 
the maneuver design are given in reference 10. 
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Instrumentation and Data Processing 


The shuttle is fully instrumented and has a number of redundant systems for mea- 
suring various vehicle states and controls. The instrument packages will he men- 
tioned specifically. First is the Aerodynamic Coefficient Tdentif ica tion Package 
( ACIP ) , an ins trumentation package specifically designed to measure rates, and accel- 
erations and control surface positions required for parameter identification. The 
ACIP data were recorded at 170 samples per second. Second is the instrumentation for 
the flight quidance and control system, the Rate Gyro Assembly, and Accelerometer 
Assembly (RGA, AA) , which is a source for acceleration and rate measurements. The 
RGA, AA data are recorded at 25 samples per second but is very noisy. The third 
source of flight measurements is the navigation instrumentation, the Inertial Mea- 
surement Unit (IMU). The IMU measurements are high fidelity but are recorded at only 
one sample per second which limits their usefulness. 

The ACIP data are the primary source for linear and angular accelerations, angu- 
lar rates, and control surface deflections. However these data were fully available 
only for flights STS-6, 7, and 8. On flights STS-1 1 and 13 the yaw rates failed to 
be recorded; an attempt to compensate for this loss was made by incorporating RGA yaw 
rate measurements. However in this study better estimates for over half the maneu- 
vers on the two flights were found using IMU rather than RGA-corrected ACIP measure- 
ments. On flight STS-1 7 a power loss resulted in a failure of any ACIP data to be 
recorded; parameter extraction was based solely on RGA, AA measurements. For all the 
flights, RCS chamber pressures were used to determine jet thrust; these measurements 
came from the vehicle operational instrumentation. 

The data considered most reliable were used to generate a best estimated trajec- 
tory (BET) for the shuttle vehicle. The data written to tapes for the parameter 
extraction consisted of only those maneuvers considered appropriate for extraction. 
The linear and angular rates and control surface deflections came from the ACIP 
instrumentation except as noted. The BET angular rates and linear accelerations at 
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the start of a maneuver were taken as initial conditions, and the rates and accelera- 
tions were integrated over time to obtain angular positions and vehicle velocities. 
The velocities were then corrected for the effect of winds, and the resulting compo- 
nents were used to calculate the vehicle total velocity, angle of attack, and angle 
of sideslip. This combined data set is recorded at 25 samples per second and com- 
prises the data contained on the tape to be processed by the parameter extraction 
software. Additional details on the instrumentation and data processing can be found 
in references 11, 12, and 13. 


Equations of Motion 

The lateral-direction equations of motion used in this study are based on per- 
turbations about trimmed flight conditions and are written relative to the body axes 
shown in figure 1 . The equations are 

• " • 

0 = fc + 8. ) + ^ cos 0 sin <j> + p sin a - r cos a (1 ) 

mv y u v 


. I XZ . 1 Y I % r xz qSb „ 

p = ■— r + — + — PP + T~ c £ 

x x x x 


( 2 ) 


. I XZ . T X - X Y I X7, qSb _ 

r = p + — pq - — qr + — C n 

Z Z Z 7j 


(3) 


<j> = p + r cos <J> tan 0 + sin <p tan 0+4) 


(4) 


where 


C = C +C 8 + C + C + C ( fir - fir ) 

Y Y Y„ Y 2V Y 2V Y ' t 

o 0 p r 6r 


+ C (6a - 6a ) + C 6RCS 

6a 6RCS 


(5) 
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( 6 ) 


( 7 ) 

The results of this study are based on maneuvers performed at velocities of 
Mach 1 and higher. For this reason the terms containing velocity are sufficiently 
small that the equations of motion are considered essentially insensitive to the 
rotary derivatives and to C. . and C •; therefore, these derivatives are fixed 

P P 

at zero throughout this study. 

Time histories of five measured quantities were fit during the estimation 
process. These are the sideslip angle (p), roll and yaw rates (p,r), lateral 
acceleration (a y ), and bank angle ($). 

Maximum Likelihood Estimation 

Stability and control derivatives were extracted using the maximum likelihood 
estimator. Among other statistical properties, the maximum I i kel ihood est imator 
is efficient and asymptotically unbiased. This estimator consists of maximizing 
the likelihood function of the measurement errors, which is the product of the 
probability density functions evaluated at each measurement time. This approach 
requires that the form of the measurement error distribution is known; it is 
normally assumed this distribution is Gaussian. 

It is assumed the actual system can be modeled as « 

• 

x(t) = F(X,U,Q,t) (8) 


C X = c x * c t e * c t W * C 1 w * v !v * c t < 6 ' • «V 

x x o x 0 x p r x 0 x 6r 1 


* C Y (6a - 6a t ) ♦ C Y 6RCS 
&a 1 6RCS 


C=C ♦ C 8 C ^ ♦ C ~ ♦ C . ♦ C (6r-6r) 

n n o V % 2V "r 2V n 0 2V n 6r * 


♦ C (6a - 6a t ) ♦ C 6RCS 
6a 6RCS 
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i 


(9) 


Y(t.) 


G(X,M,Q,t. ) + e. 


1 


where equation (B) is a vector representation of equations (1) to (4) and equation 
(9) is a vector representation of the measurements. In these equations , X is the 
state vector, U the vector of controls, Q the vector of stability and control 
derivatives, t is time, and e^ the vector of measurement noise for the measure- 
ments at time t^ • 

If it is assumed that the measurement noise is Gaussian, then the likelihood 
function (ref, 14) is 


L(Y,p) = [(2ir) 4 R]“ k/2 exp{- 1 J [Y^t.) - Y(t.)] T R _1 ^(t.) - Y( t ± ) ] } (10) 


where the subscript M denotes actual measurements and R is the measurement co- 
variance matrix. Takinq the natural logarithm of equation (10) and multiplying by -1 
yields the cost function 


J (0) = -log L(Y,p) 


I \ IVV 





N 

2 


log R + 2N log 2 tt 


(11) 




Maximization of equation (10) with respect to 0 is equivalent to minimization of 
equation (11) with respect to Q. The last term on the right is constant relative to 
p and can be neglected; if r is known, the second term can also be neqlected for 
the same reason. Minimization of the remaining term results in solving VJ * = 0 


which gives the estimates 


Q=P 




Vi 



[v 2 j(p J ]" 1 VJ(P ), 


j - 0,1,2,... 


( 12 ) 
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Since a sequence of estimates, Cb , are obtained iteratively, the process must beqin 

A 

with initial parameter estimates, Q . 

o 

If R is unknown in equation (11), direct minimization of J(Q) with respect 
to Q and R is complicated by the fact that R is an implicit function of Q. A 
simpler approach is to minimize with respect to Q and R independently. Minimiza- 
tion of equation (11) with respect to R yields 


* -W.\ (w - (VV - '"Vl 1 

1 = 1 


(13) 


The procedure used here is, first, assuming R is diaqonal with initial esti- 
mates for the diaqonal elements, iterate equation (12) several times. Then, on each 

A A 

succeeding iteration, first estimate R using the most recent value of Q in equa- 


tions (9) and (13), and then apply equation (12) once using R in J(Q) . This two 
step process is repeated each iteration to convergence. 

The computer software used to obtain the maximum likelihood estimates is MMLE3 
(ref. 14). A detailed description of the software can be found in the reference. 


Analysis and Results 

In this section the results obtained in this study are discussed. These results 
are based on extracting the stability and control derivatives from 38 maneuvers on 
the six flights. The time span for the measurements obtained during the maneuvers 
ranged from 4 to 1 5 seconds with the measurements sampled 25 times a second. 

The estimation approach taken here is based on information contained in measured 
accelerations and rates, various trajectory parameters and the measured atmosphere. 
The method of analyzing atmospheric measurements which accounts for spatial, diurnal, 
and semidiurnal corrections is described by Price (ref. 15). This atmospheric infor- 
mation is combined with onboard measurements of accelerations and rates in order to 
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construct the trajectory (ref. 16) which is used for estimating the stability and 
control derivatives. 

In the results presented, moment derivatives are relative to the fliqht center 

of qravity and were estimated with rotary derivatives fixed at zero and C y fixed 

6a 

at the data book value of 0.00042 per deqree. All mass properties and center of 
gravity information were supplied by NASA Johnson Space Center and are shown in ta- 

A 

ble 1. The weighting matrix (inverse of the measurement noise covariance matrix, R) 
was initially set to a diaqonal matrix with the values 796.3, 234.8, 4324, 237.5, and 
21820. These values correspond, respectively, to the measured variables 0, p, r, 4>, 

A 

and a y . Estimation of R using equation (13) beqan on iteration 4 for each maneu- 
ver; from 8 to 20 iterations were required for convergence. 

The extracted stability and control derivatives will be presented in figures as 
functions of Mach number. Roth flight-extracted and predicted values along with 
variations associated with the predicted values will be shown. For example, fiqure 2 
shows rolling moment due to sideslip anqle as a function of Mach number with the pre- 
dicted values (P) and variations (V) indicated by solid lines, the extracted values 

; 

by the symbol . The predicted values are based on data book values, corresponding 
to fliqht 7, which are the result of numerous preflight tests of shuttle aerodynamics 

(ref. 1). The variations reflect uncertainties in the data book values; they are 

based on differences between flight and predicted results for previously researched 
aircraft and extrapolated to the shuttle configuration. 

Lateral-Directional Moment Derivatives 

C — Extracted values of the rolling moment due to sideslip are shown in fiq- 
ure 2. Except for a few outliers, the values fall within the variations. Above Mach 
7 the flight results are slightly more positive than the predicted values, showing 
less stability than predicted. Similar results have been reported by Maine and Iliff 
(ref. 17) and Kirsten et al. (ref. 18). The estimates in the region above Mach 22 
are generally based on maneuvers having low dynamic pressure (q < 10 psf), making it 
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difficult to estimate stability and control derivatives. 'This circumstance may par- 
tially account for the estimates lying outside the variation band. 

Below Mach 7 the estimates are highly scattered. At the lowest Mach numbers, 
both aileron and rudder controls are simultaneously active. As presently configured, 
it is not possible to perform maneuvers which allow isolated control surface motions, 
thus making it difficult to accurately separate the effects of different surfaces. 
Significant differences in extracted coefficients have been noted between values when 
estimating rudder parameters versus not estimating rudder parameters for the same ma- 
neuver (ref. 4). Furthermore, as the figure shows, the uncertainty in the estimates 
grow dramatically below Mach 3. The outliers below Mach 5 occurred on fliqhts 11 and 
17. Generally, therefore, results below Mach 5 must be accepted with caution. 

Similar results were obtained with Columbia flight data (ref. 9). 

C — Results for the yawing moment due to sideslip are shown in figure 3. 

0 0 

This coefficient is similar to the rolling moment due to sideslip in that there is 
considerable scatter below Mach 7 and the estimates generally lie within the varia- 
tion band above Mach 7. This coefficient tends to be less negative than predicted 
below Mach 5 and more negative with a general downtrend above Mach 7. The large out- 
lier near Mach 1 and the outlier at Mach 14 occurred, respectively, on STS-17 and 
STS- 1 3. 

Lateral Control Derivatives 

C — Figure 4 shows the results for the rolling moment due to aileron. Below 
oa 

Mach 7, the aileron tends to be less effective than predicted; above Mach 15, aileron 
effectiveness tends to be greater than predicted. The three outliers in the lower 
left corner of the figure were extracted from STS-17 data (all measurements from the 
RGA, AA, and no PTI maneuvers). In general, aileron effectiveness tends to increase 
with increasing Mach number. 

C -- In general the coefficient of yaw due to aileron (fig. 5) tends to be 
n 6a 

less effective than predicted, although almost all of the extracted values lie within 
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the variations. The positive values below Mach 3 are highly suspect in view of the 


large uncertainty at the low Mach numbers. 

The lowered effectiveness of both aileron derivatives is consistent with the 
Columbia results in reference 9. This conclusion is especially true for both the 
Challenger and the Columbia orbiters below Mach 7. 

C — The rolling moment due to rudder is shown in fiqure 6. Almost all of 

Q 

*6r 

the estimates lie within one variation of the data book values and show this deri- 
vative to be close to what was predicted. Since most values are less than the data 
book values, there is a suggestion that the rudder may be somewhat less effective 
than predicted, especially below Mach 2.5, The negative outlier at Mach 1 was ex- 
tracted from flight 17 data, 

C — Figure 7 shows the yawing moment due to rudder. Most of the flight 
n 6r 

values lie within one variation of the data book value. However, all of the values 
also indicate the rudder to be less effective than predicted. For both rudder deriv- 
atives the Challenger results confirm the Columbia results (ref, 9) which showed 
these derivatives to be less effective than predicted. 

Side Force Derivatives 

C -- Generally, the side force derivatives are slightly more difficult to 
Y Sa 

estimate because the signal input to the estimation program has a very small signal 
to noise ratio. In addition, force signals tend to look the same regardless of 
cause, and hence, it is difficult for the program to decompose the signal into caus- 
ative components. Thus, since C is very small (0.00042) compared to other force 

6a 

derivatives, it was not possible to get a consistent estimate of this derivative with 

high confidence. Further, C appears to alias the RCS side force derivative when 

6a 

it is estimated. Therefore, for all cases presented in this report C was fixed 

6a 

at the data book value. 

C -- Side force derivative with respect to sideslip angle is shown in fig- 

Y e 

ure B. Of the ten outliers, seven are from the flights (11, 13, and 17) for which 
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AC IP information was missing. Five of the outliers occur below Mach 2 which is a 
region of great uncertainty; except for the two positive values, the remaining three 
values may be reasonable. Most values are moderately scattered within the variation 
bounds. Roth the outliers and the values within the variation bounds tend to be more 
positive than the data book values. This suggests the shuttle vehicle is less stable 
than the data book indicates, in agreement with reference 9. 

C — The side force due to rudder given in figure 9 indicates a considerable 
Y 6r 

scatter in the estimates below Mach 2 where there is great uncertainty* These re- 
sults may be indicative of the aforementioned small signal to noise ratio in the on- 
board accelerometers and the ensuing difficulty in decomposing the signal. Above 
Mach 2 the values are close to the predicted but indicate the rudder to be less 
effective than predicted. Compared to the Columbia results (ref. 9) which were 
highly scattered, the Challenger values show a definite trend. 

RCS Derivatives 

The RCS jets were treated in MMLK3 as if they were an additional aerodynamic 
control surface. The solutions were obtained throughout the speed range for side 
force, rolling-moment, and yawing-moment derivatives due to yaw jet firings. In this 
paper, yaw jet evaluation is presented as a function of Mach number on a per jet ba- 
sis. Comparisons are made to STS-7 preflight values based on known vacuum thrust 
corrected for altitude effects. Because the altitude profiles of the six flights are 
slightly different, the flight values will differ somewhat from the preflight values 
presented here. Fur thermore, the preflight values have not been corrected for flow- 
field interactions . 

C — Side force due to yaw jet firing is shown in figure 10. The differ- 

y rcs 

ences between predicted and flight values can be attributed to jet- interaction ef- 
fects consisting of flow-field interactions and vehicle impingements, in addition to 
the aforementioned altitude profile differences. The figure shows good agreement 
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between fliqht and predicted values with an indication that the yaw jets are somewhat 


more effective than predicted. 

C -- The flight values for the yawing moment due to yaw jets shown in fiq- 

n 

r cs 

ure 11 qenerally aqree well with the predicted values. Considering the sources of 
differences noted previously, the yaw jets are apparently less effective than pre- 
dicted by not more than 10 percent. The lowered effectiveness is particularly 
evident in the Mach 10 to 20 range. 

C — In the case of the rolling moment due to yaw jets shown in figure 12, 

£ rcs 

the differences between fliqht and predicted values are siqnif icantly larqer. This 
suqqests greater interaction effects than seen in the previous two derivatives. The 
greater scatter in this derivative across the Mach range indicates there is also much 
more variability in the interactions. Verification of the interactions at a few 
points using the Development Flight Instrumentation (DFI ) is given in references 4 
and 9. Thus, it appears that the lower effectiveness of this derivative can be 
largely attributed to flow-field interactions which were not originally modeled in 
the data book values. 

Overall, the RCS derivatives extracted from the Challenger flights are compara- 
ble to those obtained from the Columbia flights. That is, for both flights the same 
RCS derivatives are less/more effective over the same Mach ranges. 

CONCLUDING REMARKS 

The lateral stability and control of the shuttle orbiter Challenger has been 
analyzed over the hypersonic speed range of Mach 1 to Mach 25. Acceleration and rate 
measurements made during 38 lateral maneuvers on flights 6, 7, 8, 11, 13, and 17 were 
used in a maximum likelihood estimation program to extract the aerodynamic coeffi- 
cients. The flight-derived coefficients were compared to preflight data book values 
and previously obtained values from flights of the Columbia shuttle vehicle. 
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The extracted stability and control derivatives were usually within one varia- 


tion of the pref light values, although the scatter is generally greater below Mach 5. 
Several coefficients were found to be somewhat less effective than predicted; this is 
particularly true for the aileron derivatives below Mach 7. The yaw jet results show 
these jets to be fully effective reqardinq side force. On the other hand, the yaw 
jets appear to be only about 90 percent effective in terms of the yawing and rollinq 
moments. For the latter derivative, the lower effectiveness is apparently due to 
flow-field interactions . All of the conclusions obtained from the Challenger data 
agree with and reinforce those obtained previously from Columbia data (ref. 9). 
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TABLE 1 


ENTRY PHYSICAL CHARACTERISTICS OF SPACE SHUTTLE CHALLENGER 


Mass properties (range for six fliqhts): 
Mass , kg 


86,514 - 93,191 


Moments of inertia (range for six flights): 

l x , kg-m 2 

l y , kg-m 

I rjr f kg-m “ 

J XZ' kc *" m? 


Wing : 

Reference area, m 2 

Mean aerodynamic chord, m 
Span, m 


1,201,401 - 1,224,002 
8,904,347 - 9,435,717 
9,306,246 - 9,844,278 
176,584 - 202,778 

249.91 

12.06 

23.79 


Elevon (per side): 

Reference area, m 2 19.51 

Mean aerodynamic chord, m 2.30 


Rudder (per side panel); 

Reference area, m 2 9.30 

Mean aerodynamic chord, m 1.86 
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Figure 2.- Rolling moment due to sideslip versus 
Mach number 
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Figure 3.- Yawing moment due to sideslip versus 
Mach number 
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Figure 5.- Yawing moment clue to aileron versus 
Mach number 
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Figure 6.- Rolling moment due to rudder versus 
Mach number 
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Figure 8.- Side force due to sideslip versus 
Mach number 
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Figure 9.- Side force due to rudder versus 
Mach number 
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